A new framework for studying spherically symmetric static solutions in f(R) gravity 
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We develop a new covariant formalism to treat spherically symmetric spacetimes in metric f(R) 
theories of gravity. Using this formalism we derive the general equations for a static and spherically 
symmetric metric in a general /(iZ)-gravity. These equations are used to determine the conditions 
for which the Schwarzschild metric is the only vacuum solution with vanishing Ricci scalar. We also 
show that our general framework provides a clear way of showing that the Schwarzschild solution is 
not a unique static spherically symmetric solution, providing some incite on how the current form 
of Birkhoff 's theorem breaks down for these theories. 

PACS numbers: 



I. INTRODUCTION 

Ever since the publication of the Schwarzschild solu- 
tion almost a hundred years ago, the study of spheri- 
cally symmetric solutions has played a fundamental role 
in determining understanding of the nature of gravity 
and underlies many of the key tests of Einstein's theory 
of General Relativity. Until recently, General Relativity 
has been unquestionably the only theory able to explain 
gravity on both astrophysical and cosmological scales. 
With the advent of new high precision cosmological tests, 
capable of probing physics at very large redshifts, this 
situation has completely changed and recently the large- 
scale validity of General Relativity has begun to be ques- 
tioned. This is largely due to the fact that in order to 
fit the standard model of cosmology, which is based on 
General Relativity coupled to standard matter (baryons 
and radiation) , the introduction of two dark components 
are needed to achieve a consistent picture. Specifically, 
dark matter is needed to fit the astrophysical dynam- 
ics at galactic and cluster scales, while a new ingredient, 
dubbed dark energy, is required in order to explain the 
observed accelerated behavior of the Hubble flow. Com- 
bining the luminosity distance data of Supernovae Type 
la [l|, the Large Scale Structure @, the anisotropy of 
Cosmic Microwave Background [|| and Baryon Acoustic 
Oscillations [4] suggest that if we retain General Relativ- 
ity as the theory of the gravitational interaction, the best 
fit model is a spatially flat Universe, dominated by cold 
dark matter (CDM) and Dark Energy (DE) in the form 
of an effective cosmological constant. Although CDM 
candidates have not yet been directly detected, there 
are strong arguments that suggest that CDM has a non- 
gravitational origin The same is not true for DE. The 
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cosmological constant and coincidence problems together 
with the fact that there are no convincing DE candidates, 
seems to suggest that the Concordance model is incom- 
plete, and despite enormous effort over the past few years, 
this problem remains one of the greatest puzzles in con- 
temporary physics. One of the theoretical proposals that 
has received a considerable amount of attention recently 
is that Dark Energy has a geometrical origin. This idea 
has been driven by the fact that modifications to Gen- 
eral Relativity appear in the low energy limit of many 
fundamental schemes [f| 0] and that these modifications 
lead naturally to cosmologies which admit a Dark Energy 
like era |cHL2| without the introduction of any additional 
cosmological fields. Most of the work on this idea has 
focused on fourth order gravity, in which the standard 
Hilbert-Einstein action is modified with terms that are 
at most of order four in the metric tensor. 

The features of fourth order gravity have been ana- 
lyzed with different techniques [13j and all these studies 
suggest that these cosmologies can give rise to a phase of 
accelerated expansion, which is considered to be an im- 
portant footprint of Dark Energy. In particular, dynam- 
ical systems analysis shows that for FLRW models there 
exist classes of fourth order theories which admit a tran- 
sient decelerated expansion phase, followed by one with 
an accelerated expansion rate. The first (Friedmann-like) 
phase provides a setting during which structure forma- 
tion can take place, followed by a smooth transition to a 
DE-like era which drives the cosmological acceleration. 

More recently, the theory of linear perturbations for 
these models has been developed [11] using the 1+3 co- 
variant approach [l5| , A number of important features 
were found which allows one to differentiate the structure 
growth scenario from what occurs in General Relativity. 
Firstly, it was found that the evolution of density per- 
turbations is determined by a fourth order differential 
equation rather than a second order one. This implies 
ffialptilMsoeSotaaiaifloctt ke zappstcy. dlui@bu®tiam: aa: it ains , in 
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general, four modes rather that two and can give rise to 
a more complex evolution than the one of General Rel- 
ativity (GR). Secondly, the perturbations are found to 
depend on the scale for any equation of state for stan- 
dard matter (while in General Relativity the evolution of 
the dust perturbations are scale- invariant). This means 
that even for dust, the evolution of super-horizon and 
sub-horizon perturbations are different. Thirdly, it was 
found that the growth of large scale density fluctuations 
can occur also in backgrounds in which the expansion 
rate is increasing in time. This is in striking contrast with 
what one finds in General Relativity and would lead to 
a time- varying gravitational potential, putting tight con- 
straints on the Integrated Sachs- Wolfe effect for these 
models. 

Although these results are very encouraging, there are 
still some important open problems to be addressed. Of 
particular interest is the degree to which the physics of 
fourth order gravity is consistent with both cosmological 
and solar system scales, indeed, there has been consid- 
erable debate about the short-scale behavior of higher 
order over the past few years leading to much work on 
the Newtonian and post-Newtonian limits of these theo- 
ries [l7j . Consequently, measurements coming from weak 
field limit tests like the bending of light, the perihelion 
shift of planets and frame dragging experiments repre- 
sent critical tests for any theory of gravity. Fundamen- 
tal to confronting such tests with fourth order gravity 
is the existence of physically viable spherically symmet- 
ric solutions in these theories. The aim of this paper is 
this therefore twofold. Firstly to determine a set of gen- 
eral results for spherically symmetric spacetimes in f(R) 
gravity, and secondly to obtain a general procedure for 
generating solutions of this type. 

The present analysis is based on a powerful extension 
of the 1+3 covariant approach in which the three spatial 
degrees of freedom are further decomposed relative to a 
spatial vector [l8[ . In the case of spherical symmetry, this 
is chosen to be the radial direction. This leads to a larger 
set of covariant variables with their corresponding equa- 
tions (evolution, propagation and constraint). Further- 
more all the equations are developed in the Jordan frame 
without resorting to any conformal transformations. 

Unless otherwise specified, natural units (h = c = 
= 8nG = 1) will be used throughout this paper, Latin 
indices run from to 3. The symbol V represents the 
usual covariant derivative and d corresponds to partial 
differentiation. We use the —,+,+,+ signature and the 
Riemann tensor is defined by 

-K bed = 1 bd,c — 1 bc,d + 1 bdt ce — 1 bet de , (,-LJ 

where the T a i,d are the Christoffel symbols (i.e. symmet- 
ric in the lower indices), defined by 

= \g ae {9be,d + 9ed.b - 9bd,e) ■ (2) 

The Ricci tensor is obtained by contracting the first and 



the third indices 

Rab — g Cd Racbd ■ (3) 

The symmetrization and the antisymmetrization over the 
indexes of a tensor are defined as 

T(ab) = 7J {Tab + T ba ) , T[ ab ] = - (T ab — T ba ) . (4) 

Finally the Hilbert-Einstein action in the presence of 
matter is given by 

A=^Jd i x^g~[R + 2£ m ] . (5) 

II. GENERAL EQUATIONS FOR FOURTH 
ORDER GRAVITY 

In a completely general context, a fourth order the- 
ory of gravity is obtained by adding terms involv- 
ing f(R,R ab R ab ,R abcd R abcd ) to the standard Einstein 
Hilbcrt action. However, we know the Gauss-Bonnet 
term {Q = R 2 - AR ab R ab + R abcd R abcd ) is a total dif- 
ferential in four dimensions and hence do not affect the 
field equations. Hence we can replace all linear terms of 
RabcdR abcd with the other two. Furthermore, if the space- 
time is highly symmetric, then the variation of the term 
R a bR ab can always be rewritten in terms of the variation 
of R 2 [llj[2(|. It follows that a sufficiently general fourth- 
order Lagrangian for a highly symmetric spacetime only 
contains powers of R and we can write the action as 

A = \ J d^x^f—g [f(R) + 2C m ] , (6) 

where C m represents the matter contribution. 

Varying the action with respect to the metric gives the 
following field equations: 

f'G ab =T%+~(f-Rf)g ab +V b V a f'-g ab V c V c f, (7) 

where /' denotes the derivative of the function '/' w.r.t 
the Ricci scalar and T" b L is the matter stress energy tensor 
defined as 

T a n b = fi m u a u b + p m h ab + Cu b + ql n u a + ttS- (8) 

Here u a is the direction of a timelike observer, h ab is 
the projected metric on the 3-space perpendicular to u a . 
Also fi m , p m , q m and 7r^ denotes the standard matter 
density, pressure, heat flux and anisotropic stress respec- 
tively. Equations (0 reduce to the standard Einstein 
field equations when f(R) = R. 

III. 1+1+2 COVARIANT APPROACH 

We know that the 1+3 covariant approach, initially de- 
veloped by Ehlers and Ellis [l5[ has proven to be a very 
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useful technique in many aspects of relativistic cosmol- 
ogy. The approach has been particularly useful in obtain- 
ing a deep understanding of many aspects of relativistic 
fluid flows, whether it is applied in terms of fully non- 
linear GR effects or the gauge invariant, covariant per- 
turbation formalism. In cosmology these methods have 
been applied for example to the formalism and evolu- 
tion of density perturbations [l6| in the universe and to 
the physics of cosmic microwave background 121]. This 
approach is based on a 1+3 threading decomposition of 
the spacetime manifold w.r.t a timelike congruence as 
a splitting of spacetime onto a timelike and a orthog- 
onal three-dimensional spacelike hypersurface. All the 
essential information in the system is captured in a set 
of kinematic and dynamic 1+3 variables that have a well 
defined physical and geometrical significance. These vari- 
ables satisfy a set of evolution and constraint equations 
derived from the Bianchi and Ricci identities, forming a 
closed system of equations for a chosen equation of state 
describing matter. 

A natural extension to the 1+3 approach, optimized for 
problems which have spherical symmetry, is the 1+1+2 
formalism developed recently by Clarkson and Barrett 
fl8| . In this formalism one first proceeds to the same split 
of the 1+3 approach and then a further one that isolates 
a specific spatial direction. This allows us to derive a set 
of variables that are more advantageous to treat systems 
with one preferred direction. For example in spherically 
symmetric system the equation for the 1+1+2 variables 
are scalar equations and are much simpler than the ones 
of the 1+3 formalism which are in general tensorial. The 
1+1+2 formalism was applied to the study of linear per- 
turbations of a Schwarzschild spacetime [18j and to the 
generation of electromagnetic radiation by gravitational 
waves interacting with a strong magnetic field around a 
vibrating Schwarzschild black hole [22j |. 

In the following we give a brief review of these for- 
malisms, before applying it to the specific case of f(R)- 
gravity. 



A. Kinematics 

In (1+3) approach first we define a timelike congru- 
ence by a timelike unit vector u a . Then the spacetime 
is split in the form R ® V where R denotes the timeline 
along u a and V is the 3-space perpendicular to u a . Then 
any vector X a can be projected on the 3-space by the 
projection tensor h% = + u a u b . 

At this point, two derivatives are defined: the vector u a 
is used to define the covariant time derivative (denoted 
by a dot) for any tensor T a ' b c ..ti along the observers' 



worldlincs defined by 



rpa..b „,&X7 rpa..b 

J- c.d — U V e J c ..d 



(9) 

and the tensor h ab is used to define the fully orthogonally 
projected covariant derivative D for any tensor T 



a..b 



with total projection on all the free indices. Angle brack- 
ets to denote orthogonal projections of vectors and the 
orthogonally projected symmetric trace-free PSTF part 
of tensors 23j: 



y(a) _ n a b y b ; rp{ab) _ 



h^ r .h b K. 



-h ab h, 



cd 



(11) 

In the (1+1+2) approach we further split the 3-space V, 
by introducing the unit vector e a orthogonal to u a so that 



e a u a = , 
Then the projection tensor 

Nj EE h n b 



e a e 



g a b + u a u b ■ 



e a e a = 1. (12) 
e a e b , N a a — 2 , (13) 



projects vectors onto the 2-surfaces orthogonal to e a and 
u a , which, following [l8j], we will refer to as 'sheets'. 
Hence it is obvious that e a N ab = = u a N ab . As we 
know in (1+3) approach any second rank symmetric 4- 
tensor can be split into a scalar along u a , a 3- vector and 
a projected symmetric trace free (PSTF) 3-tensor. In 
(1+1+2) slicing, we can take this split further by split- 
ting the 3- vector and PSTF 3-tensor with respect to e a . 
Any 3- vector, ip a , can be irreducibly split into a compo- 
nent along e a and a sheet component orthogonal to 
e a i.e. 

$ a = ^ e a + ^J a , $ = ^ a e a , * a ee N ab ^ b . (14) 

A similar decomposition can be done for PSTF 3-tensor, 
ipab, which can be split into scalar (along e a ), 2-vector 
and 2-tensor part as follows: 



Ipab = 

where 



ip(ab) = * (e a e b - \ N ab^j + 2* (a e b) + f afc , (15) 



N b e c ^ bc , 



1p{ab} 



N c {a N b) d - -N ab N cd ) V> cd , (16) 



and the curly brackets denote the PSTF part of a tensor 
with respect to e a . 
We also have 

2 

h{ab} = , A <Qb ) = -e( e 6 ) = N ab - -h ab . (17) 

The sheet carries a natural 2- volume element, the alter- 
nating Levi-Civita 2-tensor: 



&ab — £abc& 



Vdab c e c u d 



(18) 



D P T 



a..b 



h a f hP c ...h b q W d h r e V r Tf-\ 



(10) 



where e abc is the 3-space permutation symbol the vol- 
ume element of the 3-space and r\ a bcd is the space-time 
permutator or the 4- volume. 

With these definitions it follows that any 1+3 quan- 
tity can be locally split in the 1+1+2 setting into only 
three types of objects: scalars, 2-vectors in the sheet, and 
PSTF 2-tensors (also defined on the sheet). 
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B. Derivatives and the kinematical variables 

Apart from the 'time' (dot) derivative, of an object 
(scalar, vector or tensor) which is the derivative along 
the timelike congruence u a , we now introduce two new 
derivatives, which e a defines, for any object ip a ...b c '" d '- 

4> a .. b c - d = e f D f ip a .. h c ~ d , (19) 
S f ^ a .. b c - d = Nj...N b g N h c ..N i d N f :i D j ^ f .. g i - j .(20) 

The hat-derivative is the derivative along the e a vector- 
field in the surfaces orthogonal to u a . The S -derivative 
is the projected derivative onto the sheet, with the pro- 
jection on every free index. We can now decompose the 
covariant derivative of e a in the direction orthogonal to 
u a into it's irreducible parts giving 

D a eb = e a a b + -<f>N ab + & ab + Cab , 



where 



Cab 



5 a e a 
1 



(21) 

(22) 
(23) 



(24) 
(25) 



^e ab 5 a e b , 
5{ a e b } ■ 

We see that for an observer that chooses e a as special 
direction in the spacetime, represents the expansion of 
the sheet, Cab is the shear of e a (i.e. the distortion of the 
sheet) and a a its acceleration. We can also interpret £ as 
the vorticity associated with e a so that it is a represen- 
tation of the "twisting" or rotation of the sheet. 

Using equations (fT4"|) and (|T5|) one can split the (1+3) 
kinematical variables and Weyl tensors as 



Ae a 
fle a - 



-A a 

n a 



(26) 
(27) 



\N ab 

t" {<.,< I, ~ \ N ab 



a ab = £ e a e b 



E r ,b — 



2E 



(a e 6) 



28 



(a^b) 



h Sab 
&ab 5 



H e a e b 



\N ab 



(28) 
(29) 

2rl {a e b) + rl ab . (30) 



where E ab and H ab are the electric and magnetic part of 
the Weyl tensor respectively. Therefore the key variables 
of the 1+1+2 formalism are 

{e,A,n,Z,£,rl,A a ,n a ,X a ,£ a ,rl a ,Zab,£ab,rlab} • 

(31) 

Similarly, we may split the anisotropic fluid variables q a 
and n ab : 



q a = Qe a + Q a , 



TTab = II 



1 



e a eb 



-N ab 



(32) 

2n (a e b) + n ab . (33) 



The full covariant derivatives of e a and u a in terms of 
these variables are given in Appendix [XJ 



IV. 1+1+2 EQUATIONS FOR LRS-II 
SPACETIMES 



Locally Rotationally Symmetric (LRS) spacetimes 
posses a continuous isotropy group at each point and 
hence a multi-transitive isometry group acting on the 
spacetime manifold [24|. These spacetimes exhibit lo- 
cally (at each point) a unique preferred spatial direction, 
covariantly defined, for example, by either vorticity vec- 
tor field or a non- vanishing non-gravitational acceleration 
of the matter fluids. The 1+1+2 formalism is therefore 
ideally suited for covariant description of these space- 
times, yielding a complete deviation in terms of invariant 
scalar quantities that have physical or direct geometri- 
cal meaning (25[. The preferred spatial direction in the 
LRS spacetimes constitutes a local axis of symmetry and 
in this case e a is just a vector pointing along the axis 
of symmetry and is thus called a 'radial' vector. Since 
LRS spacetimes are defined to be isotropic, this allows 
for the vanishing of all 1+1+2 vectors and tensors, such 
that there are no preferred directions in the sheet. Thus, 
all the non-zero 1+1+2 variables are covariantly defined 
scalars. The variables, {.4, 9, 0, £, E, fi, £ , H, fi,p, II, Q} 
fully describe LRS spacetimes and are what is solved for 
in the 1+1+2 approach. A detailed discussion of the co- 
variant approach to LRS perfect fluid space-times can be 
seen in [24j . 

A subclass of the LRS spacetimes, called LRS-II, con- 
tains all the LRS spacetimes that are rotation free. As 
consequence in LRS-II spacetimes the variables fi, £ and 
% are identically zero and the variables 

{A, 9, </>,£,£, ii,p,Il,Q} 

fully characterize the kinematics. The propagation and 
constraint equations for these variables are: 



Propagation equations 





4> =- 


1 9 






2^ + 
2 

"3^" 


£ - 








In =- 

2 





2 
1 

- -( 
3 



(34) 
(35) 

(36) 
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Evolution equations 



This gives the 3-Ricci-scalar as 



.A 



= - + 



, 1 1 



+^(^ + 3 P )-£ + in , 



+ , (37) 



(38) 



|E-e) fi + i(E-|e)n 



5 i? 



3 

4< 



A' 



(46) 



where A is the Gaussian curvature of the sheet, 2 i? a b = 
KN ao . From this equation and (jMJ) an expression for A 
is obtained in the form 25 1 



(47) 



From <|34fl39[) . the evolution and propagation equations 
of A can be determined as 



Propagation/Evolution Equations 



A 



e - e a 



A = -</>A. 



(48) 
(49) 



fi+Q = 

-p + ii = 



- e( M + P )-(^ + 2^)Q--sn{40) 



a n 



-e + s q 



i - e = - (^ - 



i 



+ 1(^ + 3?) 



Commutation relation 



ip — ip = —Aip + 



E V 



(41) 



(42) 



(43) 



From equation (|48|) . it follows that whenever the Gaus- 
sian curvature of the sheet is non-zero and constant in 
time, then the shear is always proportional to the expan- 



sion as E 



e. 



Let us now turn to the case of spherically symmetric 
static spacetimes which belong naturally to LRS class II. 
The condition of staticity implies that the dot derivatives 
of all the quantities vanish. Furthermore the expansion 
also vanishes, as a non- vanishing expansion would imply 
that the timelike congruence would contract or expand in 
time, which is not possible in a static spacetime. Hence 
we have = 0, and as discussed in the previous section 
this imply E = 0. From equation (|37p we then have 
the heat flux Q to vanish identically in these spacetimes. 
Hence the set of (1+1+2) equations which describe the 
spacetime become, 



Here the quantities \i and p are the total effective energy 
density and pressure. In context of fourth order gravity 
we would define these quantities later. Since the vor- 
ticity vanishes, the unit vector field u a is hypersurface- 
orthogonal to the spacelike 3-surfaces whose intrinsic cur- 
vature can be calculated from the Gauss equation for u a 
that is generally given as [25j |: 



^Rabcd — (Rabcd)± ~ K ac Kbd + Kb c K a d 



(44) 



where ^R a bcd is the 3-curvature tensor, _L means pro- 
jection with h a b on all indices and K ao is the extrinsic 
curvature. With the additional constraint of the vanish- 
ing of the sheet distortion £, i.e. the sheet is a genuine 

2- surface, The Gauss equation for e a together with the 

3- Ricci identities determine the 3-Ricci curvature tensor 
of the spacelike 3-surfaces orthogonal to u a to be 







P + Ii 
A 



^2 2 



1 



-u- -U-£ 

3 P 2 



- 3 A £+ l n 



(50) 
(51) 



Acj>+^(ji + 3p)-£ + ^U ,(52) 

3 
2 



A ) H- (n + p)A , (53) 
1 



{A + 4,)A+-{^ + ip) 



(54) 



V. SPHERICALLY SYMMETRIC STATIC 
SPACETIMES IN HIGHER ORDER GRAVITY 



Rab — — 









e a eb- 







-cb 2 -K 



N ab , (45) 



At this point one can re-derive these equations in the 
case of /(A)-gravity. The quantities /i, p and n are de- 
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fined, in this case, as 
1 



1 



W - /) + fx 



n 



/' V 2 

+/"A> + /"'A 2 ) (55) 

j(p m + \{f-Rf)-\f"x- 

lf»X<l>-lfi"X 2 -Af"x\ (56) 



where we have defined R = X. We will consider here the 
"external" field generated by a point-like source so that 
[i m = and p m = 0. Because of the additional degrees 
of freedom the equations (|50M54p are not closed and we 
have to add an additional equation, the Trace equation: 

Rf - 2/ = -3f"X - Zf"X<t> - 3f"X 2 - dAf'X . (58) 

Using the above equations in (|50H54p and eliminating £ , 
we get the set of four coupled first order equations gov- 
erning the spacetime in the fourth order gravity as 



+f"X (<j> + 2A) , 
/' [A + A(A + <j>) 
R = X . 



A 



1 



Rf - 77/ 



fx 



1 



Rf 



f-f"x 2 



(59) 

f'X A, (60) 
(61) 

X{4> + A)f" .(62) 



We emphasize here that the above system of equations 
are written in terms of the covariant quantities in the 
1+1+2 splitting and absolutely co-ordinate independent. 
Note that the system reduces to the second order system 
of GR in vacuum [18|], if we put f(R) = R, R = and 
X = 34] . However as in the case of the Einstein equa- 
tion, or any other fully covariant system of equations, 
the physics can be understood fully only if one chooses 
an observer. In the 1+3 approach this is done basically 
choosing a velocity field, but in the 1+1+2 framework 
this is not sufficient. One has to give also a particular 
form of ' radial co-ordinate. This in turn will define a 
specific form for the l ha£ derivative. As we will see in 
the later sections there is a natural choice for this coor- 
dinate given by the geometry of our problem and we will 
use it to find exact spherically symmetric solutions for 
some specific f(R) gravity models. 



VI. COVARIANT RESULTS FOR THE 
SPHERICALLY SYMMETRIC SYSTEM 

From the structure of (|59H62[) we can already de- 
duce some important results for spherically symmetric 



static solutions in a general f(R) gravity in absolute co- 
ordinate independent manner. These results are impor- 
tant because they can be used as guidelines to understand 
the behavior of any proposed f(R) model in this setting 
and to design new ones. 



A. Necessary condition for existence of solutions 
with vanishing Ricci scalar. 

It is evident from the equations (|59M62I) above, the 
function / must be of class C 3 at R — 0, which imply, 



|/'(0)| < +00 , |/"(0)| < +00 , |/"'(0)| < +oc 
Also we impose the conditions 

/(0) = 0,i? = 



(63) 



(64) 



Note that the condition of vanishing of the Ricci scalar 
throughout the manifold automatically implies X = 0. 
Now there are two possibilities: 

a - /'(0) 7^ : In this case we see the system re- 
duces to the following: 



A + A(A - 



-A]=0 







(65) 
(66) 



It can be easily checked that the conditions R = 0, 
/(0) = and /'(0) + 1 makes the Einstein Tensor G a b 
vanish and therefore Schwarzschild solution is the only 
spherically symmetric static solution. This then allows us 
to state a generalization of Birkhoff 's Theorem in higher 
order gravity. 

For all functions f(R) which are of class C 3 at R = 
and /(0) = while f(0) + 1 0, Schwarzschild solution 
is the only static spherically symmetric vacuum solution 
with vanishing Ricci scalar. 

It is also interesting to note that the above result is 
consistent with the conditions /' > and /" > 0, which 
guarantee the attractive nature of t he g ravitational inter- 
action and the absence of tachyons 12j . This shows that 
there may be a connection between this solution and the 
very nature of the gravitational interaction. 

The presence of this solution, can have interesting con- 
sequences on the validity of these models on the Solar 
System level. In particular if one concludes that the Sun 
behaves very close to a Schwarzschild solution, the exper- 
imental data of the solar system would help constraining 
these models. 

b. /'(0) = 0,/(0) = : In this case (15911621) are 
identically satisfied for all values of <fi and A that guar- 
antees R — and hence X = [35[ . Hence for all models 
with / (0) = 0, any solution with vanishing Ricci Scalar 
in General Relativity would be a solution to the above 
system. This is interesting as it shows that fourth order 
gravity in this context can present the same solutions of 
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GR plus additional. For example, the Reissner Nord- 
strom solution which represent the space time outside 
a spherically symmetric charged body, is a solution to 
the system (|59H62p even if no electric charge is present. 
Similarly a static spherically symmetric solution for a 
perfect fluid interior with equation of state p = (l/3)/o 
(for example Haj j-Boutros solution or the special case of 
Whittaker solution [27[) can be a solution of this system 
in the absence of any standard fluid. 

The presence of solutions of type (6) shows that when 
the conditions given in paragraph (a) are not satisfied the 
Schwarzschild solution is not a unique static spherically 
symmetric solution. Such results hint towards a disproof 
of the general Birkhoff theorem in its classical form for 
fourth order gravity. 



B. Necessary condition for existence of solutions 
with constant scalar curvature 

Solutions with constant Ricci scalar are characterized 
by the fact that R = Rq = const, and, as consequence, 
X, X = 0. Imposing these conditions on (|59H62j) and 
supposing it to be at least of class C 3 in R = Rq one 
obtains 



Jo 



A 



6 



1 



fo A + A{A + 
= -Raft + 2/o , 
where f'(Ro) = fo e ^ c - A first solution exists if 
f^O, fo^O, 2/0-^0/0 = 0. 



(67) 

(68) 
(69) 

(70) 



Instead in the case /q 7^ 0, fo = one obtains again 
the Schwarzschild solution (i?o = 0). Finally another 
solution can be achieved if 

fo = , fo = , R = Ro, X,X = 0, (71) 

is satisfied. As in the previous subsection, in this case 
also, any constant Ricci scalar solution in GR would iden- 
tically be a solution to the system. The relation (|70l) was 
already found by Barrow and Ottewill [28||_in the cosmo- 
logical context and later rediscovered in 29]. It relates 
the value of the constant Ricci scalar with the universal 



constants in the action. For example if we have the La- 
grangian as R — 2A, which is the Lagrangian for GR with 
the cosmological constant, we must have, as well known, 
the relation Rq — 4 A. 



C. The curious case of R gravity. 

As we have already explained, the condition for exis- 
tence of solutions with covariantly constant scalar curva- 
ture connects the constant curvature with the universal 
constants of the Lagrangian. However this is not the case 
for f(R) = KR? . In fact for this type of Lagrangian the 
third condition of (|70l) is identically satisfied. This means 
that we can have a constant curvature solution for any 
value of the curvature. Thus for R 2 gravity, the 'cosmo- 
logicaF constant term in a Schwarzschild-dS/ AdS space- 
time becomes a local constant of integration just like the 
mass. Hence in this theory we can have two distant stars 
behaving like two different Schwarzschild-dS/ AdS object 
with different values of the constant. Unfortunately this 
case is rather pathological since it corresponds to the 
case in which the trace of the field equations in vacuum, 
3D/' + f'R — 2/ = is satisfied identically for constant 
Ricci scalar, whereas usually it may be satisfied for spe- 
cial values of R. In any case this model is ruled out by 
solar system experiments (see [Hj]). 



VII. CHOOSING A CO-ORDINATE SYSTEM 
AND RELATION BETWEEN THE C O VARIANT 
VARIABLES AND THE METRIC 

The most natural way to choose the proper radial co- 
ordinate in spherically symmetric static spacetimes, is to 
make the Gaussian curvature 'K' of the spherical sheets 
to be proportional to the inverse square of the radius. In 
that case, this co-ordinate V becomes the area radius of 
the sheets. This gives a geometrical definition to the 'hat' 
derivative. As we have already seen, K = —<pK, therefore 
the most natural way to define the hat derivative of any 
scalar M would be 



(72) 



With this choice, the system of equations (|59H62j) be- 
comes 
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f ( \r4& + h 2 -A^) = LRf If 



dr 



,., ,1 dA 



2 
A 2 



2A)f"X, 



" 4 "' = \ f ~\ Rf '~ Af " X 



1 ,dR 
J 2 y dr 



X 



-Rf' + -f-f"'X 2 -(Xcf> + XA)f". 



(73) 

(74) 

(75) 
(76) 



r 



In the r co-ordinate above the most general spherically 
symmetric static metric is 

ds 2 = -A{r)dt 2 + B(r)dr 2 + r 2 {d6 2 + sin 2 9d(f> 2 ) (77) 

Now, from the properties of the four-velocity u a and the 
radial vector e a i.e. u a u a = — 1 and e a e a = 1, we find 
that 



11 i = V / iW,e r = v / BM, (78) 

also, from the definitions of different covariant scalars, 
we get 



A = -u b u a V b e a 



1 dA 
2A~d^ 



D 



(79) 



provided that n — 1, 2, > 3. However since R — is by 
itself a differential constraint involving and A, hence 
any <fi and A that ensures a zero Ricci scalar would solve 
the system. As we know the following solution 



r 



2m . m 
, ■A- - T 



1 - 



2m 



(85) 



with equations (|T9"|) and (|50|) gives the usual 
Schwarzschild metric in (t, r, 0, 0) co-ordinates that 
has a zero Ricci scalar, hence the above solution is the 
solution of the system. 



N\V b e a = -y/B 
r 



(80) 



R can be found in the usual way as a contraction of 
the Riemann tensor and X is derived from it as in (|75|) . 
Thus we see any solution to the equations (I73H76[) . would 
uniquely determine the metric for the spacetime. 



B. A solution with constant non-zero Ricci scalar 

As described before if we substitute X = 0, R = Rq ^ 
in the above system of equations then a solution is 
possible if and only if n = 2. In that case the solutions 
of the other two functions are 



VIII. AN EXAMPLE: SOME EXACT 
SOLUTIONS FOR R n GRAVITY 

In this section we present, as an example, a few exact 
solutions for R n gravity, in absence of standard matter. 
Specializing the choice of f(R) = R n , equations (|73H76[) 
becomes 



2 dr 



1 dA _,_]. 
-nrd)——R 

2 dr 



A 



1 



4>R r ' 



-R n 



1 dR 
—r 
2 

dX 



—rm 

2 V dr 



-r(fm{n - ^)~^R 



n-2 



2 J 3n 
+{n - l)R n - 2 X {<j) + 2A) (81) 

1 — 2T7 

-(A + (f>) AR"- 1 + — — R n 
on 

-{n-l)R n - 2 XA, (82) 

X , (83) 
2-n 



-R n - X 



A) 



-n{n - l)(n - 2)R n ^X 2 (84) 
A. Schwarzschild solution 

Substituting R = 0, dR/dr = in the above set of 
equations, we see that the equations are satisfied trivially 



r V t 3 



A 



m + R r 2 



2m R Q 2 

1 + — r 

r 3 



(86) 
(87) 



This is the usual Schwarzschild-dS / AdS solution depend- 
ing on the sign of Rq. 



C. A solution with non-constant Ricci scalar 
vanishing at infinity 

To find more non-trivial solutions of the above system 
of equations, let us use a Schwarzschild like ansatz, 



cp = VCir" + C 2 r? , R = C 3 /r~< ( 7 > 0), 



such that the Ricci scalar vanishes at infinity. We use 
these ansatz in the system of equations and then alge- 
braically solve for the powers and coefficients such that 
the system is identically satisfied. With this choice we 
get the following solution: 
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A 



-C (5 - 4n) r 4 " + (An 2 - 6n + 2)r _1 

2(2 - n) 

/ 

(l + 2n- 2n 2 ) (7 - 10n + 4n 2 ) 



^(l + 2n-2n 2 ) (7-10n + 4n 2 ) (l + CY~ 7 ""i 4 " ) ^ 



V 



i? = 



X = 



{2-nf (l + Cr 

6n(n — 1) 



(2n (n - 1) - 1) r 2 ' 

12n (n - 1) [ (l + 2n - 2n 2 ) (7 - lOn + 4n 2 



(2n(n-l)-l)r3 I (2 _ n)2 ^ + Cr -^±^) 
Now solving for the metric coefficients, we get 



A(r) =r {2n - 2) i^T + 



(2n-l) (7 



(2 - n) 2 



(2-n)* 



' B(r) (7-10n + 4n 2 )(l + 2n-2n 2 )) 



C 



(7-10re + 4ti' ! ) 



(99) 

(90) 
(91) 
(92) 



(93) 



This solution was originally found by Clifton 30]. The 
solution reduces to Schwarzschild for n = 1 and valid for 
n < (1 + v3)/2 beyond which the metric has unphysical 
signature. However for n G (1,(1 + \/3)/2) the Ricci 
scalar is negative and hence the action is only real- valued 
if n is an even rational number. That is in it's lowest form 
the numerator of the fraction is even. This problem can 
be also be resolved by assuming the absolute value of 
the Ricci scalar in the action. However in that case the 
Schwarzschild limit at n = 1 is not possible as \R\ is not 
differentiable at R — and hence doesn't belong to the 
class C 3 functions. It is also interesting to note that in 
spite of the Ricci scalar vanishing at infinity this solution 
is not asymptotically flat. 

IX. CONCLUSION 

In this paper we have analyzed static spherically sym- 
metric metrics within the /(i?)-gravity framework. Using 
the 1+1+2 formalism we were able to derive equations 
describing these metrics for a general form of the func- 
tion / and a pointlike source. These equations have been 
used to obtain a set if general conditions for the existence 
of certain types of static spherically symmetric solutions. 
It is important to note that our system of equations are 
much simpler than what one obtains when writing the 
Einstein Field Equations in terms of metric components 
and therefore it is much easier to find new solutions and 
general covariant results. In particular the results when 
/(0) = and /'(0) = are much more transparent using 
this approach. 

Our results show that the presence of solutions with 



constant Ricci Scalar is influenced by the properties of 
the derivatives of the function / up to the third order. 
This implies that two f(R) models are indistinguishable 
in this framework if differences only arise after the third 
derivative of / at the given value of the Ricci scalar. 
Also one can probe in general that a form of the Birkhoff 
theorem exists for f(R) gravity only if /(0) = and 
/'(0) 7^ 0, while in general there is more than a static and 
spherically symmetric solution for the field equations. 

It is also interesting to note that the conditions /' > 0, 
/" > 0, which guarantee the attractive nature of the 
gravitational interaction and the absence of tachyons [l2[ 
are consistent with the the form of Birkhoff's theorem 
stated here. This suggests a possible link between these 
conditions and the Birkhoff theorem, something which 
definitely deserves further study. 

In order to extract observable results the 1+1+2 equa- 
tions need to be further specialized choosing a specific 
form of the radial coordinate. This is equivalent to the 
choice of an observer in the 1+1+2 formalism, which, un- 
like the 1+3 case, requires not only the specification of a 
velocity field but also a specific spatial direction. The re- 
quirement that the Gauss curvature has an inverse square 
dependence offers a natural choice for this coordinate. 

Once this is done given any /(i?)-theory of gravity (and 
sufficient ingenuity) one can derive static and spherically 
symmetric solution(s) for this theory. We have used a 
i?™-gravity as an example to derive some exact solutions. 
As expected from our general considerations, since for 
this class of models /'(0),/(0) = 0, the system ([5911521) 
does not admit a unique solution and Birkhoff's theorem 
is violated. 

It is worth to stress, however, that such considerations 
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are limited to the case of pointlike sources. It is known 
that the situation can be really different int he case of 
extended ones [33| . Such issues will be treated elsewhere. 

As a final comment we would like to point out that if 
one admits in the fact that in this framework the Birkhoff 
theorem is violated, any Newtonian limit of a background 
solution will give, in principle, a different results [26| . 
This means that there is no way to calculate the physical 
Newtonian potential without knowing the exact back- 
ground which characterizes the entire Universe. This is 
not surprising because in these theories the relation be- 
tween local physics and the rest of the Universe is much 



tighter that in General Relativity due to their relation 
with Mach's principle. 

Appendix A: Useful relations 

In this appendix we give some useful relation needed 
for the calculations performed in the text. 

The full covariant derivatives of e a and u a are 



Y„< i, -- -Au a u b - u a a b + ( E + -6 j e a u b + (E a - e ac Q, c ) u b + e a a b + ^<t>N ab + £e ab + Cab , (Al) 



V a ti h 



-u a {Ae b + A b ) + e a e b ( -6 + E 



+ e a (E fc + e hc O c ) + (E a 



i e b - 



-&£ab 



(A2) 



The covariant time derivative of e a is given by, 

e a — Au a + a a , where A — e a u a , (A3) 

and a a is the component lying in the sheet. 

The new variables a a , 0, £, Cab and a a are fundamen- 
tal objects of the spacetime and their dynamics give us 
information about the spacetime geometry. The spatial 
covariant derivative of a scalar ^ is defined as 



and 



(A4) 



while for any vector ^ a orthogonal to both u a and e a 
(i.e. \I/ a lies in the sheet), the various parts of its spatial 
derivative may be decomposed as follows (Note that a bar 
on a particular index indicates that the vector or tensor 
lies in the sheet.): 



D a * b = -e a e b ^ c a c - e b 
+ e a i>- b + S a ^ b . 



-6^ a + {£,£ac + Cac) * ' 



(A5) 



Similarly, for a tensor ty a b (where ^ a b = ^{ab}) '■ 

Sa^bc 

-4>^ c)a + V c ) d (& ad + Cad) 



D * bc = -2e a e [b ^ c)d a d + e a ^ bc 
-2e 1 



(b 

For the Levi Civita 2-tensor we have: 

£a&e 6 = = £( ab ) , 

^abc — ^a^bc ^b^-ac T" 6 c £ ab 

£abS Cd = N a C N b d - N a d N b C , 

£a C £bc = N ab , £ ab £ab = 2 . 



•(A6) 



(A7) 
(A8) 
(A9) 
(A10) 



£ a b — ~2u[ a £ b ] c A c + 2e[ a e b ] c a c 
iab = 2e[ a e b j c a c , 
S c £ab = . 



For the projection tensor we have 



(All) 



Nab 


= 2u (a ?i b) 


N ab 


= -2e (a a 


ScNab 


= . 



(A12) 
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